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SECTION 1 


INTRODUCTION 

This final report documents the work completed by the 
University of Dayton Research Institute under a grant with the 
Dryden Research Facility of NASA Ames Research Center. The 
purpose of the grant was to further develop and evaluate the 
FASTEX flutter analysis computer code. Modifications and 
improvements to the code have been documented and are presented 
herein . 

1.1 OVERVIEW 

After a discussion of the general background of flutter 
analysis and the FASTEX code, Sections 2, 3, and 4 of this report 
document the relevant equations of motion, modal interpolation 
procedures, and control system considerations, respectively. 
This information is followed in Section 5 with a summary of the 
software developments. Additional information documenting input 
instructions, procedures, and details of the plate spline 
algorithm are found in the appendices. 

The system of codes delivered to NASA is referred to as 
UDFASTEX . A magnetic tape, listings and sample runs are 
submitted under separated cover. These codes will allow NASA to 
make more efficient use of computer resources in the analysis of 
advanced aerospace vehicles. Key features of the codes are their 
improved modularity, more extensive internal documentation, and 
ease of interfacing new aerodynamic methods, and reduced user 
intervention in the interpolation of modal data. The UDFASTEX 
system operates on the VAX 11/700 series computers and makes use 
of user-friendly interactive menus written for operation under 
the VMS operating system. 
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1.2 BACKGROUND 


In the development of advanced aerospace vehicles, it is 
essential to perform sufficient engineering calculations to 
assure that the vehicle will operate safely within its flight 
envelope. This is especially true for the major airframe 
components such as the wings, canards, and empennage. 
Calculations need to be performed to verify that these key 
airframe components have adequate strength and stiffness to 
withstand both the static and dynamic loadings. 

Lifting surfaces must be designed to be aeroelast ically 
stable and yet light enough to achieve overall performance 
requirements. It was with that objective in mind that the FASTOP 
computer code was originally developed in 1975 [1]. Both 
strength and flutter calculations are performed in FASTOP with 
the goal being to resize the structure so that near minimum 
weight can be achieved without violation of the constraints on 
stresses, deflections, and flutter speeds. 

The FASTOP code was further improved since 1 975 by the 
University of Dayton [2-4] and by Grumman Aerospace [5] by 
improving the optimization procedure and its structural analysis 
capability. However, by far the greatest use of the code has 
been to calculate flutter and divergence speeds. Its 
aeroe last ic ity capability has a number of state of the art 
features which make it a highly versatile flutter analysis tool. 

In 1 980 the Air Force Flight Dynamics Laboratory awarded a 
contract to the University of Dayton to streamline the FASTOP 
flutter analysis capability. The primary objective was to update 
the code so that it would operate only as an analysis tool 
without the optimization capability. The code was reorganized 
for faster execution times on the CYBER-175 system at Wright- 
Patterson Air Force Base. Additional features such as a more 
general modal interpolation procedure were also included. The 
resulting code, FASTEX [6], has since been extensively used in 
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both government and industry for the analysis of advanced flight 
vehicles such as the X-29A. 

Although the code is now operational the VAX 11/780 
minicomputer, further improvements were necessary to remove 
problem size limitations for operation on the CYBER. Modal 
interpolation between a structure grid and an aerodynamic grid 
have been improved to take advantage of the VAX virtual memory. 

1.3 OBJECTIVES 

The objectives of this grant have been to increase the 
problem size capacity of FASTEX, reduce run times by modification 
of the modal interpolation procedure, and to add new user 
features. Interfaces have been provided to aid in the inclusion 
of alternate aerodynamic and flutter eigenvalue calculations. 
Other user features have been added to increase the ease of 
program use. 

1 . 4 SUMMARY 

All modifications to the FASTEX code are operable on the VAX 
11/700 series computers under the VMS operating system. The 
program can now be dimensioned to desired problem size and 
compiled. Aerodynamic and flutter solution methods are now 
operable as separate packages. Plots can be made of the flutter 
velocity, display and frequency data. A preliminary capability 
has also been developed to plot contours of unsteady pressure 
amplitude and phase. 
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SECTION 2 


EQUATIONS OF MOTION 

In this section the aeroelastic equations of motion are 
developed in modal coordinates. This is presented as basis for 
reference for the code modifications presented in Section 5. 

Following the development of Meirovitch [7], we apply 
Hamilton's principle in the form 


(2-1 ) 



where 

6W = -6V+5W^c = virtual work 
6T - variation of kinetic energy 

6V = variation of potential energy 

SW„„ = virtual work due to nonconservative forces. 
NC 

The modal assumptions 


N 

(2-2) w - l Mx,y ) q. (t ) 

i-1 


N 

(2-3) P = I P, (x,y ) q. (t) 

i-1 1 1 


are made where 


w = w(x,y,t) = lateral deflection at point x,y at time t 
p - net upward air load at point x,y at time t 
4 k - assumed deflection mode 
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= pressure contribution due to <Jk mode 
q^ » generalized coordinate. 

This leads to the kinetic energy expression 


(2-4) T = 1 { q } T [ M] {q } 

where the elements of the [ M ] matrix are 

(2-5) M. . = J1 m(x , y ) 4>.(x,y) <J> .(x,y) dS 

1J S 1 J i-1 , N; j-1 , N 

with 

m = mass per unit area of the lifting surface 
S = surface area. 

Also, the potential or elastic strain energy takes the form 
(2-6) V = 1 { q } T [K] {q } 


where the elements of the [K] matrix are 

(2-7) K. . = \\ k(x,y) 4> i ( x , y ) <Mx,y) dS 

1 J q 1 J 


and 

k(x,y) = stiffness influence function. 

The final term in the formulation requires the variation of 
the nonconservative work term which has the form 

(2-8) 6W nc = { 6q } T [P] {q} 

where the terms in the pressure matrix are 

(2-9) P.j - Jj f g 4> i (x ,y ) p j (x ,y ) dS 
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Taking the variation of the kinetic and strain energies and 
using Equations 2-8 and 2-9 in Equation 2-1 leads to the equation 
of motion 

(2-10) CM] {q } + [K] {q} - [P] {q} = {0} 

Introduction of the assumption of harmonic motion and 
hysteretic damping gives 

(2-11) -oi 2 [M]{q} + ( 1 +ig) [K] {q } = [P]{q} 

where the [P] matrix is determined using one of the conventional 
unsteady aerodynamic programs such as doublet lattice. The [P] 
matrix is expressed as 

(2-12) [P] = ^pV 2 [Q] 

and substitution into (2-11) yields 

(2-13) A[K] {q } - (CM] + [Q*]) (q) 

where 


A = ( 1 + i g ) / io 2 

[Q*] - c [Q ] 

1 ,,2 . 2 

c = g P V /to 

Note that the reduced frequency is 
k = <ub/V 
so [Q*] becomes 

(2-14) [Q*] * (c/1 /2 p V 2 ) [ P] 
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Since [Q] is the generalized force per unit dynamic pressure 

(2-15) [Q] = d [Q*] 

2 2 

where d = 1 /c = 2k /pb . The aerodynamic forces output by the 
UDFASTEX routine called MATDAT. 

The flutter eigenvalue problem solved is 
(2-16) [ A] {q } - A[ B] {q } 

where 

[A] » [M] + [Q*] 

and 

CB] = [K] 

Since [Q*] is a function of reduced frequency, k, and Mach 
number, M, the [A] matrix is also 

[A] = [ A( k , M ) ] . 
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SECTION 3 


INTERPOLATION PROCEDURES 


Several situations arise in a flutter analysis where it is 
necessary to perform interpolations on given or computed data. 
These interpolations generally involve data of either one or two 
variables. The discussion in this section is directed toward the 
following three common interpolation problems: 

• Evaluation of aerodynamic influence coefficient matrices 
at intermediate values of reduced frequencies for a fixed 
Mach number (i.e., one-dimensional interpolation). 

• Evaluation of aerodynamic influence coefficient matrices 
at intermediate pairs of reduced frequency and Mach 
number . 

• Calculation of modal deflections and slopes at the 
aerodynamic grid points in terms of data based on the 
structural model. 

The first interpolation problem alone is an essential feature of 
the p-k method. Interpolation of modal data is required in 
nearly all problems with the exception of piston theory 
calculations. In the piston theory case, aero and structural 
grid points can be easily made identical. The use of doublet- 
lattice or Mach box methods requires a specialized grid which 
rarely is the same as the structures grid. 

In the following subsections, spline and polynomial 
interpolation procedures are discussed. It is generally 
considered that splines yield less oscillatory interpolated 
values than polynomial approximations based on least-squared 
fits. Polynomial approximations are discussed since many of the 
air loads programs (see References 8 and 9, for example) require 
as input the coefficients to polynomials which approximate the 
mode shapes. 
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3.1 ONE-DIMENSIONAL SPLINE 


A one-dimensional cubic spline is a piecewise polynomial 
defined on an interval (x^.x^). Functional values are assumed to 
be known at a given set of points x.j ,X£ . . • • » x ^ which are called 
knots. The polynomial is so constructed that its functional 
values equal the given values at the knots and furthermore its 
first and second derivatives are continuous at these points. The 
work by deBoor [10] provides much useful information regarding 
both the theoretical and practical aspects of spline interpo- 
lation . 

To evaluate aerodynamic matrices at intermediate values of 
reduced frequencies, use can be easily made of IMSL routines 
[11]. For a given real or imaginary term in an aerodynamic 
matrix, its value is known for a given number of reduced 
frequencies 

jkg* • • • *k • It is desired to determine the value of this term 

for some value of reduced frequency, k, which lies between k^ and 

k . This is accomplished by two IMSL calls. The first is to the 
n 

routine ICSCCU which obtains the values of the coefficients to 
the spline polynomial. The second call is to ICSEVU which uses 
the coefficients to compute the given term at the selected value 
of k. Loops must be set up to evaluate all terms, both real and 
imaginary, in the aerodynamic matrix. Numerical experimentation 
using the IMSL routines, together with Reference 8 doublet 
lattice program, has shown that quite accurate results can be 
obtained for a given set of widely spaced k values. 

The above procedure can also be used to evaluate aerodynamic 

matrices when both reduced frequency and Mach number are treated 

as variables. Assume that k^,...,k are the given reduced 

frequencies and that M.,...,M are the Mach numbers. Thus the 

1 m 

aerodynamics are known for all pairs (k.,M.) for i-1 n and 

j-1,...,m. To find the aerodynamics at any (k,M) pair which lies 
within the (k^.Mj) grid the following procedure is suggested for 
any given real or imaginary term in the aerodynamic matrix: 
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1. Set M=M 1 and find the value of the term at k=k using the 
above IMSL routines. Save this value. 

2. Repeat Step 1 for M=M 0 M . 

2 m 

3. Using the saved values from Steps 1 and 2, find the 
value of the term at M = M using the IMSL routines, 
treating M as the variable instead of k. 

4. Repeat Steps 1, 2, and 3 for all terms in the matrix, 
both real and imaginary parts. 

By the judicious use of interpolation on the aerodynamic 
matrices, experience has shown that computing costs for a flutter 
analysis can be reduced by as much as a factor of 10. This is 
especially true for complex configurations when more than 10 
modes are used in the flutter analysis. 

3.2 TWO-DIMENSIONAL SPLINE 

The IMSL library also contains a routine, IBCCU, for 
interpolating data of two variables. Based on experience with 
interpolation of modal data given an (x,y) grid, special purpose 
procedures should be used rather than routines like IBCCU. In 
the following paragraphs a special purpose "plate spline" is 
discussed and accuracy improvement features are suggested. 

The interpolation problem is thus stated: 


Gi ven 


and 


u J ' y j' w j ) 


j*1 , . . . , N 


(x j ,y j ) . j =i , . . . , M 

fit a mathematical function W-W(x,y) through the points Wj and 
compute 
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w ( x j ,y j ) - w. 

^WCx^.yj) = W xj (optional) 

-r-~W ( x . , y . ) = W . (optional) 
oy J J A J 


The data can be coordinates and either given modal 

deflections or given modal slopes in the streamwise direction. 
It is recommended that if modal slopes are available in the 
structures grid, this should be used rather than the optional W x j 
calculations. Appendix C outlines the basic bivariate spline 
procedure based on the work of Harder and Desmarais [12]. 
Additional features that are not in the IMSL bivariate method 
are 

1. Translation of given data to a centroidal origin, 


2. Rotation of data to principal axes, 

3. Scaling of x-y data by the reciprocal of their maximum 
values , 

4. Linpack [13] is used to solve the equation. 

Use of this procedure gives quite accurate normalwash 
information for input to the aerodynamic routines. Whenever 
possible, it is recommended that structural grid data be used 
which extends as nearly as possible to the perimeter of the 
aerodynamic planform. The plate spline can extrapolate beyond 
the structures grid; however, the streamwise slopes may not be 
accurately represented. 


3.3 POLYNOMIAL MODE REPRESENTATIONS 

Programs such as the do ubl et- - 1 at t i ce [8] require that 
vibration mode amplitudes be input in the following 
nondimensional form 
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(3-D 


w^ (x ,y) 


n 

l 

i = 1 


m 

l 

j - 1 


,< k > 
1 J 


<i> 


i“1 v j-1 


where the k superscript denotes the mode number, b is the 

reference length used in the reduced frequency definition, 

( k ) 

a. . are coefficients to be determined, and n and m are integers 

X J 

which determine the degree of the polynomial in x and y, 
respectively. For a given set of r discrete deflections in the 
k fcl1 mode (<J>j k ^,...,<f>^ k ^) at the known points 
( (x^ ,y 1 (x ,y ) ) , Equation 3-1 can be put in the matrix form 


(3-2) 


U (k) } - [B] {a (k) } 


where [B] is a rectangular matrix of size r x p and p=nm. The 
coefficients of Equation 3-1 have been arranged in a column 
vector and the [B] matrix contains terms involving powers and 
products of the known point coordinates. 

In the usual case when r>p, Equation 3-2 is an over- 
determined system. In this case the standard least-squares error 
approach can be used to determine the coefficients. The IMSL 
routine for this problem is LLSQF. 

When the 1 eas t- squares polynomial approximation is used, 
great care should be taken to use as low an order polynomial as 
possible to represent the modes. The possible oscillatory 
behavior of high order polynomials can cause highly inaccurate 
values of streamwise slope which will lead to unrealistic 
pressure distributions. 
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SECTION 4 


CONTROL SYSTEM CONSIDERATIONS 

The use of active control systems to increase the 
operational range of flight vehicles which are subject to 
aeroelastic instabilities has been an emerging technology since 
the early 1970's. The signals from sensors can be used to 
command the movement of aerodynamic controls to prevent flutter 
[14]. The use of velocity-root locus plots in conjunction with 
the p-k method is a feasible approach. Both rigid body and 
flexible mode variations need to be included in the active 
control method. The work of Nissim [15-16] and Rudisill [17] is 
important from the standpoint of setting up an optimal control 
law which relates the system response to the modal deflections. 

The approach used in NASTRAN is one approach to general 
control system problems. The viewpoint is that the control 
system is accessory to the structure. The vector of structural 
degrees of freedom is augmented by the addition of extra 
variables which are particular to the control systems (e.g., 
voltages, position variables, etc.). It is assumed that these 
control system variables, denoted by u , are related to N 
structural degrees of freedom or control variables, u^ , by the 
following expression 


b Q +b 1 p+b 2 p‘ 


(i) 2 


+ a. p ) u , 


where the a's and b's are user specified constants and p is the 
first order time derivative operator. 

Much of the work on active flutter control has utilized the 
motion of flaps, tabs, and other leading and trailing edge 
control devices. It is feasible also to utilize aerodynamic 
interference effects due to adjacent surfaces. By input of 
appropriate torsional root motion to the forward surface, for 
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example, a flutter condition on the aft surface may be 
controlled . 
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SECTION 5 

SOFTWARE DEVELOPMENT 


This section describes the software modifications that were 
performed on a NASA-developed version of FASTEX. The NASA code 
used for this work is a modified form of the Reference 6 codes 
which were developed for a CYBER computer system. The primary 
NASA modifications to FASTEX dealt with file structures, program 
interfaces, and related coding changes required for operation 
under the VMS operating system on VAX 11/700 series computers. 

The system of codes described herein are referred to collectively 
as UDFASTEX . 

5.1 OVERALL APPROACH 

The overall approach taken to improving the FASTEX code was 
to reorganize the code into major program modules, provide for an 
improved modal interpolation methodology, and develop a 
convenient user interface. A single procedure manages the files 
and provides a menu-based interface for activating all the 
program options. 

Table 5.1 shows the UDFASTEX main menu. The various options 
are described in the subsections below. Appendix B documents the 
procedure file which generates this menu. When options are 
selected from the menu, further information is given concerning 
files which are needed to execute the various options. 

Figure 5.1 shows the primary program modules: 

UDFASTEX1 Collection of frequently called subroutines 

and the main calling program. This module 
is used when Option 1 is selected from the 
menu . 

UDFASTEX2 Same as above, except that it is designed 

for use with Option 2. 

POOL ET AL Routines for basic flutter input. 
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Table 5.1 


Option Selection Menu for Overall UDFASTEX System 


UDFASTEX MAIN MENU 


OPTION DESCRIPTION 

1 Run UDFASTEX with internal mode interpolation 

2 Run UDFASTEX with external mode interpolation 

3 Compute aero box data or run plate spline 

4 Run an alternate aero or eigen solution 


5 

Plot flutter V-g or 

V-f 

data 



6 

Change the problem 

si ze 

capacity of 

UDFASTEX 

7 

Compile and link an 

new 

version 

of 

UDFASTEX 


for OPTION 1 use 





8 

Compile and link an 

new 

version 

of 

UDFASTEX 


for OPTION 2 use 





9 

Create a print file 

of 

selected 

UDFASTEX output 


99 Quit 
ENTER OPTION 
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RODDEN ET_AL Subsonic aerodynamic routines. 

MIDI_ET_AL Reference 6 FASTEX modal interpolation 

rout ines . 


NEW_MIDI_ET_AL Routine which reads the downwash output 

file from the spline MODE_INTERPOLATION 
program. 


PGM_ET_AL 

SOLFLT_ET_AL 

K_S0LUTI0N 
PK SOLUTION 


Main calling routine for a supersonic 
aerodynamic code. 

Main calling routine for flutter 
eigensolution. 

k-method flutter solution routines, 
p-k method flutter solution routines. 


Table 5.2 lists the subroutines and functions included in 
these modules. The NEW_MID_ET_AL module is not shown since it 
calls no subroutines. The program modules are files which are 
compiled and linked to form the UDFASTEX1 and UDFASTEX2 
executable files. 


The calling sequence of all the subroutines 
codes shown in Table 5.2 is given in Figures 5.2 


in the UDFASTEX 
through 5.8. 
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r~ POOL ET AL 


\— RODDEN ET AL MIDI ET AL* 


UDFASTEX1 


r— PGM ET AL 


•— SOLFLT ET AL 



K_SOLUT ION 
PK SOLUTION 


*NEW MIDI ET AL for UDFASTEX2 System. 


Figure 5.1. Primary Code Modules in UDFASTEX. 
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TABLE 5.2 

SUBROUTINES AND FUNCTIONS IN UDFASTEX CODE MODULES 


Module 


Subrout ines 


Functions 


UDFASTEX1 and 
UDFASTEX2 

FOP, AFAM, AORDER ANDOR , CDABSF , 

CLUSAL , CLUSTL, COMSCA1 , CD AT , CDET , COMSCA , DCMPLF , 

CONV , DVALUE, HELGX , DSC APR , DSQRTF , F , RDM, ROUND 

ORIENT, QFLIN.RNRW 

POOL_ET_AL 

POOL, FLINFO, DENSITY 

RODDEN_ET_AL 

RODDEN, PARTI ,BIDI, MERGE, 
GLOBAL, GENQ,PRT2,QUAS, 
FUTSOL , GENF.QINTP1 , ATAN3 » 
BEIN, AUGW, INCRO , SN PDF , 
SPLIT3.TRIDI.TKER, KERNEL, 
IDF1 , IDF2 

PGM_ET_AL 

PGM 

SOLFLT_ET_AL 

SOLFLT, PRPLT, PICTUR, 
SCLMAX , SCLINC 
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TABLE 5.2 (continued) 


Module 

Subroutines 

Functions 

MIDI_ET_AL 

MIDI , MODAL, FORM, HELGA, 
HELP ,MOVIS 

- 

K_SOLUTION 

E IGM, SRE VNC , FLSL , BUCK , 

GGCHK, TRFR, FLASH, VALCOM, 
VALROM, COMVEC , CQR , CLR , CLINEQ , 
MATDAT 


PK_SOLUION 

FLOP , RTIN , ZANLYN , FRORD, 
VECP, ADI V, GCVEC.GR VEC, 
ASSESS, UERTST, GENEIG, 
JORCOM 

- 
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MAIN-FOP-AF AM 


POOL - @ 

— RODDEN - ( 5 ) 


— PGM^ 


— SOLFLT - (c) 

— PRPLT 1 


Read Input 


Subsonic Aerodynamics 


Supersonic Aerodynamics 


Flutter Solution 


Print Plots 


1 


Inactive 


in Current UDFASTEX . 


Figure 5.2. UDFASTEX Calling Sequence: MAIN. 
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"'ah of the basic flutter data (excluding aerodynamic and modal) 
are input in FLINFO. 

2 

New routine to compute air density for given altitude. 


Figure 5.3. UDFASTEX Calling Sequence: POOL. 
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Generalized aerodynamic forces are printed by GENF. 

Figure 5.4. UDFASTEX Calling Sequence: RODDEN. 
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1 

Generalized aerodynamic forces are printed by QFLIN 

2 

p-k flutter calculations 
k flutter calculations 


Figure 5.5. UDFASTEX Calling Sequence: SOLFLT. 
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1 Surface modes interpolated 

2 

Body modes interpolated 
■3 

Downwash file assembled for surfaces and bodies 
4 

Called only for UDFASTEX1 . Separate MODE_INTERPOLATION program 
used for UDFASTEX2. 


Figure 5.6. UDFASTEX Calling Sequence: MIDI, BIDI, MERGE. 
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RTIN 



ZANLYN 

FRORD ■ 
AORDER 

VECP ■ 

QFLIN - 
ADIV - 


GCVEC - 

ORIENT 

GRVEC - 
CLUTSL 


■ASSESS 

UERTST 

AORDER 

QFLIN - 

GENFIG 

RNRW 

HELGX 

CDAT 

CLUSAL 

COMSCA1 

CLUSAL 

CDAT 

COMSCA1 

CLUTSL 



HELGX 

JORCOM 

ORIENT 


Figure 5.7. UDFASTEX Calling Sequence: 


FLOP. 
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I- CONV 


h-SREVNC 


® 



I- BUCK 


GGCHK 
TRFR — 
RNRW 


—FLASH 

— C0MSCA1 

- VALCOM — 

— CQR 

- VALROM — 

— CLR 

- COMVEC 


- CLINEQ — 

— CDABSF 

- AORDER 



Creates output file with all matrices required to perform 
external flutter ei gensolution . 

2 

Prints k-flutter solution results. 


Figure 5.8. UDFASTEX Calling Sequence: EIGM. 


5.2 INCREASED PROBLEM SIZE CAPACITY 

When Option 6 is selected from the Main Menu, a Fortran 
code, SIZE, is executed which prepares a batch edit file for 
modifying the UDFASTEX "Include Files." These files are 
designated by the file type, .$NS, and they contain Key Words 
which are replaced with numerical values which determine the 
dimensioning of the program. When the batch editing job is 
completed, the Include Files have the file type, .INS. 

Table 5.3 lists the Key Words, their description, and 
default values. The SIZE program is interactive and displays the 
description and default values shown in the table. The last 
column of Table 5.3 refers to the data item and page number in 
the original FASTOP manual [1] which is still the basis for most 
of the UDFASTEX input. 

The UDFASTEX labeled common blocks are shown in Table 5.4. 

A majority of these contain the Key Word parameters which are 
replaced with numerical values. The Include Files contain 
dimensions statements and other variables which are set by the 
SIZE code. Table 5.5 is a listing of all the Include File names. 
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Table 5.3 


Key Words and Default Values for UDFASTEX Dimensions 


Key Word 

Description 

Default 

Item/Page 

$MOD 

Modes, LC( 2 ) 

20 

4/218 

$SUF 

Surfaces, LC(3) 

30 

4/218 

$REV 

Reduced Velocities, 
LC(4) 

30 

4/218 

$DEN 

Air Densities, 
LC(5) 

10 

4/219 

$MEC 

Mode Elimination 
Cycles, LC ( 25 ) 

25 

4/221 

$SVC 

Stiffness Variation 
Cycles, LC(26) 

20 

4/221 

$PBA 

Max p-k Velocities 

20 

19/227 

$PAN 

Total Panels, NP 

50 

3R/239 

$BOD 

Total Bodies, NB 

20 

3R/240 

$S0P 

Problem Size, NCORE 

400 

3R/240 

$CON 

Total Control Surfs 

5 

1 6R/246 

$M0 L 

Modal Lines, NGP 

12 

18R/247 

$PBA 

Body Points, NGP 

20 

31R/252 
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Table 5.4 


Labled Common Blocks in UDFASTEX 


BAL 

COLMEM 

DETAIL 

LABELS 

BANDW 

COMA 

DLM 

MATRIX 

BODY 

COMRWP 

DSIOI 

MESAG 

CALCP 

CONSTS 

DSI02 

MODD 

CBYTES 

CORE 

DSRN 

MODE 

CDATER 

CPLOTS 

ELMNT 

MOD V 

CEEL02 

CPLOTF 

FACE 

MODV 

CFILES 

CSETUP 

FILE 

NTPS 

CFMTA 

C SHI FT 

FITR 

PI GW 

CFMTBO 

CTABLE 

FLEXT 

PLACES 

CFMTB 

CTAPES 

FLUT 

PLAYFF 

CFMTCO 

CTFH 

FLUTAN 

PLUG 

CHEAD 

CTITLE 

FLUTB 

PRPL 

CHSP 

CTMH 

FLUTC 

REPORT 

CIDIV 

CTSH 

FLUTQ 

SIZES 

CLIST 

CTSHF 

FLUTV 

VARBLS 

CLUEF 

CTSHFO 

FSIOI 

WAYTS 

CLUEM 

CTSHV 

FSI02 

XYZ 

CLUEV 

CVIBRA 

JUNK 

YZY 

CLUFO 

DCOM1 

KLUES 


CNTRL 

DCOM2 

KLUFF 


COLDES 

DCOM3 

KMP 
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Table 5.5 


Include Files Used in UDFASTEX 


ADIV 

AFAM 

ANDOR 

AORDER 

ASSESS 

ATAN3 

AUGW 

BEIN 

BIDI 

BUCK 

CDABSF 

CDAT 

CDET 

CLINEQ 

CLR 

CLUSAL 

CLUTSL 

COMSCA 

COMSCA1 

COMVEC 

CONV 

CQR 

DCMPLF 

DSCAPR 

DSQRTF 

DVALUE 

EIGM 

F 

FLASH 

FLINFO 

FLOP 

FLSL 

FOP 

FORM 

FRORD 

FUTSOL 

GCVEC 

GENEIG 

GENF 

GENQ 

GGCHK 

GLOBAL 

GRVEC 

HELGA 

HELGX 

HELP 

IDF1 

IDF2 

INCRO 

JORCOM 

KERNEL 

MATDAT 

MERGE 

MIDI 

MODAL 

MOVIS 

NEW MIDI 

ORIENT 

PARTI 

PGM 

PICTUR 

POOL 

PRPLT 

PRT2 

QFLIN 

QINTP1 

QUAS 

RDM 

RNRW 

RODDEN 

ROUND 

RTIN 

SCLINC 

SCLMAX 

SNPDF 

SOLFLT 

SPLIT3 

SRENVC 

TKER 

TRFR 

TRIDI 

UDFASTEX1 

UDFASTEX2 

UERTST 

VALCOM 

VALROM 

V ECP 

ZANLYN 
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5.3 CALCULATION OF AERODYNAMIC PANEL DATA 


Option 3 on the Main Menu has two functions. The first is 
the calculation of aerodynamic panel and body data. The second 
is the execution of the plate spline modal interpolation. This 
subsection deals with the first Option 3 function which must be 
completed before the interpolatioh program is executed. The next 
subsection deals with the modal interpolation part of Option 6. 

The calculation of the aerodynamic panel and body data is 
accomplished by the program AER0_B0X. The panel data consists of 
the coordinates of the sending and receiving points plus the 
length and width of each panel. See Reference 8 for a discussion 
of the theory behind the aerodynamic paneling process. Reference 
18 contains the documentation for the doublet lattice computer 
code, H7WC. This code is essentially the same as the one 
currently in UDFASTEX and the previous versions of FASTEX and 
FASTOP. Figure 5.9 shows a sketch from Reference 18 and also the 
variable names for the key output of AER0_B0X. 

The input to AERO_BOX is the usual UDFASTEX_INPUT.DAT file. 
The output files created are 

AERO_BOX_OUTPUT.DAT 

AERO_BOX_XYZ_SENDING. DAT 

AERO_BOX_XYZ_RECEIVING.DAT 

AERO_BOX_LENGTH_WIDTH.DAT 

AERO_BOX_BODY.DAT 

BEIN_INPUT.DAT 

These contents of these files are described in Table 5.6. 

The AERO_BOX program contains the UDFASTEX code as the main 
program. Calls to the flutter solution routine, SOLFLT, have 
been removed. The only code from RODDEN_ET_AL that is used is 
PARTI and its associated subroutines. Subroutines RODDEN, PARTI, 
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The AERO_BOX program contains the UDFASTEX code as the main 
program. Calls to the flutter solution routine, SOLFLT, have 
been removed. The only code from RODDEN_ET_AL that is used is 
PARTI and its associated subroutines. Subroutines RODDEN, PARTI, 
and POOL are renamed AER0_B0X_R0DDEN , AER0_B0X_PART1 , and 
AERO_BOX_POOL respectively. The POOL_ET_AL routines have been 
modified so that the downwash file F0R049.DAT is not created. 

The MODE_INTERPOLATION program theory is documented in 
Appendix C. The main routine contains an overall loop on surface 
data and an inner loop on modes. The surfaces are interpolated 
first, followed by a call to BIDI_MAIN which performs the body 
mode interpolation. 

The interpolation requires as input all of the above 
AER0_B0X files except AERO_BOX_OUTPUT.DAT. The user must also 
prepare files as follows: 

MODE_POINTS_XYZ.DAT 
MODE_SHAPES_DEF . DAT 
BODY_MODES.DAT 

These files contain the given modal data and mode shape 
deflections for the surfaces and fuselage. The formats are 
included in the internal program documentation of 

MODE_INTERPOLATION. Appropriate modifications can easily be made 
to MODE_INTERPOLATION to accommodate input generated by a finite 
element code. Some sorting of the data will be required in order 
to eliminate all but the deflection components normal to the 
surface . 

Output files from MODE_INERPOLATION are FOR041 , 42 , . . . , 
FOROXX, F0R027, F0R034 where XX - 40 + number of surfaces. A 
T0TALS.DAT file is also created which contains problem size data. 
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The MODE_INERPOLATION program is intended to be run before 
execution of UDFASTEX2. All files are managed by the UD.COM 
procedure. The user is presented menus and directories 
indicating which files are needed for MO DE_I INTERPOLATION and 
UDFASTEX2 . 

5.4 OTHER PROGRAM OPTIONS 

Other program options provide for generation of flutter 
plots (Option 5), independent flutter eigensolution (Option 4), 
and compiling and linking of the various UDFASTEX versions of the 
code (Options 7 and 8). An Option 9 has also been included to 
selectively print portions of the UDFASTEX output file. This has 
not been fully implemented; however, as user requirements are 
defined for particular problem output requirements the option can 
be made fully active by studying the comments in the UD.COM file. 

The flutter solution requires the output file of the MATDAT 
subroutine and is fully compatible with the UDFASTEX file 
structures. The solution method is based on the LZ algorithm and 
solves the generalized complex eigenvalue problem without matrix 
inversion . 
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SECTION 6 
CONCLUSIONS 


A modified version of FASTEX has been developed and 
documented. The key developments are dynamic dimensioning of 
arrays and an improved modal interpolation methods for the 
lifting surfaces. 

Additional work is need to fully checkout the modal 
interpolation for the fuselage. Several users have expressed 
concern about the accuracy of the interpolated results. In 
further checking of the problem, it appears that the original 
version of FASTOP have have not been accurately coded. 

Additional work is needed to verify this. As it currently 
exists, the UDFASTEX code is fully operational for models with 
multiple lifting surfaces. It can also be used for unsymmetrical 
configurations by following the input requirements in Appendix A 
herein . 

Due to the complexity of the code and its file structures, 
further development of FASTEX is questionable. It is important 
to recoginize that UDFASTEX is a fully developed and operational 
code except as noted above. In view of new demands imposed by 
upcoming aerospace system, it appears highly desirable to 
initiate planning of a new generation of flutter analysis tools 
which will take advantage of new and efficient computer hardware. 

New code developments should pay particular attention to 
data structures and modern structured programming. Standard 
mathematical software has developed to a very high degree since 
the original FASTOP code development was planned in the early 
1970's. Utilization of this software should be planned in the 
next generation codes. 

Advances in computational fluid dynamics should also be 
reviewed. Panel methods are highly efficient; however, flutter 
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behavior in the transonic regime cannot be accurately predicted 
with current flutter analysis tools. New CFD approaches as well 
as improved panel methods need to be developed and added to a 
structured flutter analysis code. 

Alternatives to conventional eigensolut ions should be 
developed. Nonlinear effects can be handled with time 
integration of the basic equations. This can also be used in the 
next generation of codes in which active controls are included. 
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APPENDIX A 


INPUT DATA INSTRUCTIONS 


FOR UDFASTEX 


NOTE: Input changes fran previous NASA Dryden versions 
are shown in bold letters 


1.1 TITLE - 1:6 (6 lines of Title cards) 

Format (FREE) 

TITLE data are not included in UDFASTE3CL or UDFASTEX2 data 

2.1 (IC(I), 1=1,40) (Parameter inputs) 

Format (1015) 

Note: Integer parameters defined as below 

IC (1) = Defines solution type 

= 1, K type flutter analysis 
=-l, P-K type flutter analysis 
= 2, divergence analysis 
= 0, pressure calculations only 

IC(2) = Defines number of vibration modes to be used in analysis 

Maximum of 20 

IC(3) = Number of lifting surfaces (max. 30 for DIM, Doublet 

Lattice Method, 5 for MBM, Machbox Method) 

IC(4) = Number of reduced velocities (Ki) 

= 6 for P-K solution 

= 30 or less for K- solution on pressure calculation 
= 1 for divergence analysis 

IC(5) = Number of air densities (If IC(1) = 0, set IC(5) = 0) 

Maximum of 10 

IC (6) = Aerodynamic forces print option 

= 1, print data 
= 0, no print-out is required 

IC(7) = Aerodynamic pressures print option 

= 1, print data 
= 0, no print-out is required 

IC(8) = Lift and moment coefficients print option 

= 1, print data 
= 0, no print-out is required 

IC(9) = 1, reads QBAR data (frequency independent addition to 

aerodynamic matrix) 

= 0, no addition 

IC(10) = Print option for interpolated generalized aerodynamic 
forces for K- method 
= 1, print data 
= 0, no print-out is required 
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IC (11) = Index for non-zero mode numbers for normalization of 
flutter determinant. Any non-zero frequency is 
acceptable, suggest IC(ll) = 1. 

IC (12) = Index flutter determinant formulation 

= 1, for non-zero frequencies <D = K (M + A„) ' 

= 0 f in presence of zero frequencies <D = (M + A^,) K' 

IC(13) = Index for interpolation of aerodynamic forces 
= 1, perform interpolation 
= 0, to compute forces directly at each Ki 
(If IC(1) = -1, set IC (13) = 1) 

(If IC (1) = 1, set IC(13) = 0 or 1, as desired) 

(If IC(1) = 0 or 2, set IC(13) = 0) 

IC(14) = 1, option for CAICCMP plots 
= 0, no plots 

IC(15) = Index for velocity scale in flutter solution 
= 1, to use true air speed 
= 0, adopt equivalent air speed 

IC(16) = Option for addition of structural damping 

= 1, indicates addition of same damping value for all 
modes 

=-1, indicates addition of differential damping value 
for each mode 

= 0, absence of added damping 

IC(17) = Print option to display number of iteration for each 
root in P-K solution 
= 1, print data 

= 0, no print-out is required (also for IC (1) =-l) 

IC (18) = Option for root extrapolation in P-K solution 
= 1, to use two previous root values 
= 0, to use previous root value (also for IC(1)=-1) 

IC(19) = Index for ordering of roots after P-K solution 
= 1, to perform ordering 

= 0, no ordering required (also for IC(1)=-1) 

IC (20) = Print option for root iteration in flutter solution 
= 1, print iteration 
= 0, no print-out is required 

IC (21) = Index to define type of aerodynamics program 
= 1, Dm option 
= 2, MBM option 
= 3, Kernel function approach 
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IC (22) 

IC (23) 
IC (24) 
IC (25) 
IC (26) 

IC (27) 
IC (28) 

IC (29) 
IC (30) 
IC (31) 

IC (32) 
IC (33) 
IC (34) 
IC (35) 


Index for generation and storage of A^, matrix 
coefficients 

0, to compute and save 

1, assumes existence of such computed data 
-1, indicates use of MBM 

1, print modal input data (STARS/ solids data) 

0, no print-out is required 

1, print interpolated modal data (GRIDCHG output) 

0, no print-out is required 

; Number of modal elimination cycles 
Maximum of 25 

Option for stiffness variation 

1, to vary stiffness 

0, no variation required 

: Mode number for stiffness variation 
> 0, if IC (26) =0 

Option of eigenvector display 

1, print data 

0, no print-out is required 

‘ 1, display vectors 
= 0, no display 

= 1, print flutter determinant matrix for K- solution 
! 0, no display 

= Option to implanent changes in generalized masses 
and modal frequencies 
: 1 / to perform changes 
: 0, no change 

Option to revise generalized stiffness matrix 

1, to input changes 
0, no change 

= 1, for steady state analysis 
= 0, for oscillatory analysis 
(For IC (1)=2, IC (33) =1) 

= 1, option for input of scaling factor to computed 
aerodynamic force 
= 0, no factors 

= Print option for MBM 
= 1, print downwash in diaphragm region 
= 0, no print-out is required 

(If IC (21) =2 or IC (22) =1, set IC(35)=0) 
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LC (36) = Flutter redesign option 
= 1, perforin redesign 
= 0, no redesign 

LC ( 3 7 ) = Print option for DLM 

= 1, print geometric data associated with elements 
= 0, no print-out is required 

LC(38) = MATDAT output file 
= 10 

LC ( 3 9 ) = Option for future use 
= 0 

LC (40) = Density-altitude option 

= 0 for usual input of density ratios, RHOR. 

= 1 for input of altitudes (feet) in place 
of RHOR input. 

3.1 IN (Notes IN = 1 for STARS mode input by tape) 

Format (15) 

IN - 2 for input of nodal data in the input file is usable in 

UDFASTEX1 runs. The mode shape is ignored in the UDFASTEX2 run 
since that program takes input directly from the output files 
of the MOD AL__I NTERPOLAT I ON program. Include all the following 
four data items when IN 3 2 and use is made of UDFASTEX1. 

3.1.1 NO = NUMBER OF DEGREES OF FREEDOM PER MODE. FOR UDFASTEX2 RUNS 

SET TO ZERO AND OMIT INPUT OF ITEM 3.1.2. 

3.1.2 QZ (K, I) 

Modal deformations for the I-th mode and K-th degree of freedom. 

Format 7E10.0 for K = 1 to NC. Repeat for each mode, 1=1 to 
LC(2) . OMIT IF NC=0. 

3.1.3 NCARD 

Number of data lines containing all non-zero elements of the 
generalized mass matrix. Format 115. 

3.1.4 I, J, NW(I,J) 

Row, column, generalized mass matrix entry (lb.). 

Format 3 (2I5,E10.0) • Number of lines of data is NCARD. 

3.1.5 OMG(I) 

Modal frequencies (Hz.). Format 7E10.3. Number of lines of data 
is (LC (2) -1) /7+1. 

3.2 BR, FMACH (Note: BR = Reference semi-chord length, inches, 

FMACH = Freestream Mach number) 

Format 2F10.0) 

3.3 K ( I) , 1=1, LC (4 ) 

Format (7F10.0) 

Description: To read reduced velocities for K-method or 

ASE open-loop analysis 
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NV, VI, DV 
Format (15, 2F10.0) 

Note: This data is required only for P-K Solution 

NV = Number of velocities for initial analysis (maximum = 20) 
VI = Lowest velocity in knots (suggest at least 200) 

DV = Velocity increment (suggest less than 250) 
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TOL, <FVBO(I) , 1=1, 6' 

Format (7F10.0) 

Note: TOL = Tolerance for interpolation fit (usually = 10 ) 

RVBO(I) = Six reference reduced velocity values for inter- 
polation; the first and last values are to be 
within range of values specified in 3.3, for K- 
solution. For P-K method, RVBO(l) .LE.l. 69*12* 
VMIN/9BR*WMAX and RVBO(6) .GE.1.69*12*VMAX/(BR*WMIN) 


VMIN = VI knots 

VMAX = VI + (NV-1) *DV knots 

WMAX, 

WMIN = Maximum and minimum modal frequencies (rad/sec) 


3.6 MADD, I ADD, MSYM 
Format (315) 

Note: Required only of IC(31)=1 

MADD = Number of generalized mass changes 

IADD = Number of changes in modal frequencies 

MSYM = 0, if changes in generalized mass matrix are symmetric 
= 1, if such changes are not symmetric 


3.6.1 I, J, WW(I,J) 

Format (215, F10.0) 


Description: To input modification in generalized mass matrix 

Note: I,J = Row and column number of altered generalized mass 

WW(I,J) = Altered mass value (If MSYM = 0, enter only upper 
triangular masses) 

3.6.1 is repeated for K=l, MADD 


3.6.2 I, CMC (I) 

Format (15, F10.0) 

Description: Input for frequency modifications 

Note: I = Mode number to be changed 

CMG(I) = New frequency value (HZ) 

3.6.2 is repeated for K=l, IADD 


3.6.3 (MAX, GMIN, VMAX, FMAX 
Format (4F10.0) 

Note: (MAX = Minimum value of damping scale for V-g, V-f plots 
GMIN = Minimum value of damping scale for V-g, V-f plots 
VMAX = Maximum value of velocity scale in knots 
FMAX = Maximum value of frequency scale in HZ 
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3.6.4 RHOR(I) , 1=1, LC(5) 

FORMAT (7F10.0) 

NOTE: RHOR(I) * denoting ratio with reference to sea level for 

flutter and divergence analysis, a maximum of 10 

HOTE: RHOR(I) = altitudes in feet if LC(40) = 1. Conversion 

back to density ratios is Bade by call to 
subroutine DENSITY. 

3.7 NOTIR, (NINZ ( J, I) , J=l, NOTIR), 1=1, LC(25) 

Format (1015) 

Description: To read modal elimination data set, LC(25) times 

Note: NOTIR = Number of modes to be eliminated 

NINZ (J, I) = Mode number of modes to be eliminated 

3.8 FL, ACAP 
Format (2F10.0) 

Note: FL = Reference chord length ig inches 

ACAP = Reference area in (inch) 

3.9 NDELT, NP, NB, NCORE, N3 , N4 , N7 
Format (715) 

Description: Data on doublet lattice aerodynamic paneling 

Note: NDELT = 1, aerodynamics are symmetrical about Y=0 

= -1, aerodynamics are antisymmetrical about Y=0 
= 0, no symmetry along Y=0 (single surface) 

NP = total number of 'panels' on all lifting 'surfaces' 

and all optional interacting bodies, a maximum of 50. 
Each surface is divided into major trapezoidal sub- 
divisions called panels, based upon geometrical 
discontinuities. The parallel edges are parallel to 
the free stream. 

NB = number of bodies that aerodynamically interact with 
the surfaces (maximum = 20) 

NCORE = problem size, N * M 

N is the total number of aerodynamic elements 
(maximum = 400) 

M is the number of modes 

N3 =1, display pressure influence coefficients 

= 0, no display 

N4 * 1, display influence coefficients relating downwash 

on lifting surfaces to body element pressures 

N7 = 1, calculate pressures and generalized forces 

* 0, cease computation after influence coefficients 

are determined 
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4.1 XOD(I), YO(I), ZO(I) , GOIAS(I) 

Format (4F10.0) 

Description: To translate panels back to original location, perturbed 

originally in GRIDCHG program. Such coordinates are in the 
Global (Aircraft) system indicating position of origin of the 
ICS for each panel. 

Note: GG4AS ( I) = panel dihedral or V.T. rotation 

= 0 for interference panel 

4.1.1 XI, X2, X3, X4 , Yl, H2 

Format (6F10 . ) 

Description: The coordinates (XI, Yl) and (X2, Yl) are for the 

leading and trailing inboard corners of the panels. 

The coordinates (X3, Y2) and (X4, Y2) are the leading 
and trailing outboard corners in ICS. 

4.1.2 Zl, Z2, NS, NZ, COEFF 

Format (2F10.0, IX, 213, 3X, 1F10.0) 

Description: Defines boundaries of 'elements' in the panel. The 

panel is divided into a number of smaller trapezoids, 
called elanents, by lines of constant percent panel 
panel chord and of constant percent panel span. 

Note: Zl = Z-coordinate of inboard edge of the panel in inches 

Z2 = Z-coordinate of outboard edge of the panel in inches 
NS = number of element boundaries in the spanwise direction 
(maximum of 50) 

= 2, for each body interference panel 
NZ = number of elanent boundaries in the chordwise direction 
(maximum of 50) 

COEFF = entered as zero 

4.1.3 (TH(J) , J=1,NC) 

Format (6F10.0) 

Description: Chordwise elanent boundaries for the panel in fraction 

of chord (TH(1)=0.0, TH(NZ)=1.0) 

4.1.4 (TAU(J) , J=1,NS) 

Format (6F10.0) 

Description: Spanwise elanent boundaries for the panel in fraction 

of span (TAU(1)=0.0, TAU(NS)=1.0) 

Note: The above data for section 4 is to be r epeated NP times in the 

sequence as below: 

1. Vertical panels or plane of symmetry (y=0) 

2. Panels on other surfaces 

3. Body interference panels 
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4.2 XBO(I), YBO(I), ZBO(I) , 1=1, NB 
Format (3F10.0) 

Description: Input of reference coordinates, vertically vibrating 

bodies before laterally vibrating ones X, Y and Z 
global reference coordinates of the body 


4.2.1 


ZSC, YSC , NF, NZ, NY, COEFF, MRK(I,1) , MRK(I,2) 
Format (2F10.0, IX, 312, 3X, 1F10.0, 213) 


Note: ZSC 

YCC 
NF 
NZ 

NY 


COEFF 

MRK(I,1) 

MFK(I,1) 


local Z coordinate of the body axis 

local Y coordinate of the body axis 

number of body element boundaries along its axis 

1, indices body vibrating in Z-direction 

0, not vibrating in Z-direction 

1, indicates body vibrating in Y-direction 

0, indicates body not vibrating in Y-direction 
(NZ = NY) 

0.0 

index of the first elsnent on the first interference 
panel 

index of the last elsnent on the last interference 
panel 


4.2.2 (F(J) , J=1,NF) 

Format (6F10.0) 


Description: X-coordinates of the divisions of the body starting 

with body nose and proceeding aft in inches (in local 
coordinates) 


4.2.3 (RAD(J), J=l, NF) 
Format (6F10.0) 


Description: Radii of body elements at the end points of division 

in inches 
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4.3 NSTRIP, NPRl, JSPEES, NSV, NBV, NYAW 
Format (615) 

Description: Antisymmetric aerodynamics data calculations 

Note: NSTRIP = number of chordwise strips of panel elements on all 

panels 

= 1, for I£ (8) =1, used for print outs of lift and 
moment coefficients for the strips, suggest 
NSTRIP=1 

NPRl = 1, option to print pressures in QUAS or FUTSOL 
= 0, no print out is required 
JSPECS = 1, antisymmetrical aerodynamics about Z=0 
(biplane or jet effect) 

=-l , synmetrical aerodynamics about Z=0 (ground effect) 
= 0, no symmetry about plane Z=0 
NSV = number of strips lying on all vertical panels lying 
on the symmetric plane Y=0 

NBV = number of elements on all vertical panels lying 
on the plane Y=0 

NYAW =0, if NDELT=1 (symmetric about Y=0) 

= 1, if NDELT=-1 (antisymmetric about Y=0) 

= 0 or 1, if NDELT=0 

4.3.1 LIM(J,1) , LIM(J,2) , LIM(J,3) J=l, NSTRIP 
Format (313) 

Note: For NSTRIP 1, use a blank card 

LIM(J,1) = index of first element on each chordwise strip 
LIM(J,2) = index of last elanent on each chordwise strip 
LIM(J,3) = 0 

5.1 KSURF, NBOXS, NCS 
Format (115, 215) 

Description: Primary surface data associated with modal interpolation 

Note: KSURF = T this surface has one or more control surfaces with 

forward hinge lines 

= F this surface has no control surfaces 
NBOXS = number of elsnents in this surface including control 
surfaces 

NCS = number of control surfaces on primary surface 
(maximum of 5) 

Note: Repeat this data item IC(3) times (nuober of surfaces) before 
going on to 5.1.1 
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' IMPORTANT CHANGE \ 

OMIT ALL THE REMAINING ITEMS EJCEPT 5.4 IF USING UDFASTEX2 VERSION. 
MODE INTERPOLATION IS PERFOIMD EXTERNALLY USING A SURFACE SPLINE. 
CONTROL SURFACES ARE TREATS) AS SEPARATE AERODYNAMIC SURFACES. 

' IMPORTANT CHANGE \ 


5.1.1 NLINES, NELAXS, NICH, NISP 

Format (415) 

Note: NLINES = number of lines on this primary surface along which 

nodal data are input (maximum of 20) 

NELAXS * 1, translation and pitch rotation are prescribed at 
each input point 

= 0, only translation is prescribed 

NICH = control word option for the type of extrapolation done 
in the chordwise direction in interpolating modal data 
in the aerodynamic grid 

= 0, linear 

= 1, quadratic 

* 2, cubic 

NISP = control word option for the type of extrapolation done 
in the spanwise direction in interpolating modal -data 
to the aerodynamic grid 

* 0, linear 

= 1, quadratic 

= 2, cubic 

5.1.2 NGP(I) , XTEPMl (I) , YTEBMl(I), XTEPM2(I), YTERM2 (I) 

Format (15, 4F10.0) 

Note: NGP(I) * number of points on the Ith line of primary surface 

at which the nodal data are specified 
(maximum number is 12) 

XTEEMl(I) = x-coordinate of the inboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

YTERM1 (I) = y-coordinate of the inboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

XTERM2 (I) = x-coordinate of the outboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

YTERM2 (I) = y-coordinate of the outboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 
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5.1.3 (YGP(J,I), J=1,NGP (I) ) 

Format (8F10.0) 

Note: YGP (J, I) spanwise coordinates of the points along the Ith 

line at which input modal data are given in inches 
in ICS. Start with the most inboard point and 
proceed outboard. 

Note: Repeat data input in 5.1.2 and 5.1.3 for 1=1, NLINES 

5.1.4 DIST 

FORMAT (E10.0) 

Description: An arbitrary chordwise distance for a primary surface 

from the given line to a reference line on which modal 
displacements are calculated in inches (input if 
NELAXS=1) 

5.2 (XI (J), Yl (J) , X2 (J) , Y2(J), J=l, NOS) 

Format (4F10.0) 

Description: hinge line data for control surfaces. To enter data 

only if KSURF=T 

Note: XI (J) = x-coordinate of inboard terminus of the Jth control 

surface leading edge in inches in ICS 
Yl ( J) = y=coordinate of inboard terminus of the Jth control 
surface leading edge in inches in ICS 
X2(J) = x-coordinate of the outboard terminus of the Jth control 
surface leading edge in inches in ICS 
Y2 ( J) = y-coordinate of the outboard terminus of the Jth control 
surface leading edge in inches in ICS 

5.2.1 NLINES, NELAXS, NICH, NISP 
Format (415) 

Note: NLINES = number of lines on this control surface along which 

modal data are input (maximum of 20) 

NELAXS = 1, translation and pitch rotation are prescribed at 
each input point 

= 0, only translation is prescribed 
NICH = control word option for the type of extrapolation done 
in the chordwise direction in interpolating modal data 
in the aerodynamic grid 
= 0, linear 
= 1, quadratic 
= 2, cubic 

NISP = control word option for the type of extrapolation done 
in the spanwise direction in interpolating modal data 
to the aerodynamic grid 
= 0, linear 
= 1, quadratic 
= 2, cubic 
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5.2.2 NGP(I), XTEFM1 (I) , YTEEM1 (I) , XTEBM2(I), YTERM2(I) 

Format (15, 4F10.0) 

Note: NGP(I) = number of points on the Ith line of control surface 

at which the nodal data are specified 
(maximum number is 12) 

XTERMl(I) = x-coordinate of the inboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

YTERMl(I) = y-coordinate of the inboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

XTEIM2(I) = x-coordinate of the outboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

YTERM2 (I) = y-coordinate of the outboard terminus of the Ith 
line for the primary surface, in inches, defined 
in ICS 

5.2.3 (YGP (J, I) , J=l, NGP(I)) 

Format (8F10.0) 

Note: YGP(J, I) spanwise coordinates of the points along the Ith 

line at which input modal data are given in inches 
in ICS. Start with the most inboard point and 
proceed outboard. 

Note: Repeat data input in 5.2.2 and 5.2.3 for 1=1, NLINES 

5.2.4 DIST 

Format (E10.0) 

Description: An arbitrary chordwise distance for a control surface 

fran the given line to a reference line on which modal 
displacements are calculated in inches (input if 
NELAXS=1) 

5.3.1 INGP, NSTRIP, IPANE1 

Format (315) 

Description: Body surface data associated with modal interpolation 

Note: NGP = number of points on the Jth body axis at which modal 

data are prescribed (maximum of 20) 

NSTRIP = number of interference panels (or strips) associated 
with the Jth body 

IPANEL = index of the first such interference panel associated 
with the Jth body 
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5.3.2 (XGP(J), J=l, NGP) 

Format (6F10.0) 

Note: XGP(J) = streamwise coordinates of each point at which modal 

data are prescribed in inches in ICS 

Note: To repeat 5.3.1 and 5.3.2 for each body, i.e. , NB times 

5.4 KLUdB 

Format (15) 

Note: KLUGLB = 1, print global geometry 

= 0, no print out is required 
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APPENDIX B 


PROCEDURE FILE DOCUMENTATION 


PROCEDURE FILE DOCUMENTATION 


The purpose of this appendix is to document the directory 
and file structures of the main UD.COM procedure file discussed 
in Section 5. The directory structure is shown in Figure B.l. A 
description of the directories is given in Table B.l. Notes in 
this table refer to Table 5.1 which lists the UDFASTEX menu 
options . 

A listing of the UD.COM procedure is also included in this 
appendix . 


B-l 




* [DRYDEN. RFT] directory name used on University of Dayton 
VAX 11/780 


Figure B.l. Directory Structure of UDFASTEX Codes. 
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TABLE B . 1 

UDFASTEX DIRECTORIES 


Directory 

MAIN 

BIGGER 

INCLUDE 


TEST 


REVISED 


Descri ption 


Main directory for all programs, data, and 
procedure files developed. Directory 
[DRYDEN.RFT] name used on University of 
Dayton Research Institute VAX 11/780. 

Development directory used for the dynamic 
dimensioning of UDFASTEX. 

Contains all UDFASTEX Fortran codes, data, 
and procedure files for Options 1, 2, 7, 8, 
and 9 shown on Table 5.1. The Fortran codes 
contain "include" statements of the form 
Name . INS 

where "Name" is the subroutine or function 
name for the respective routines. 

Contains all UDFASTEX Fortran codes, data, 
and procedures for Option 6 shown in Table 
5.1. Files of the form 
Name . $NS 

are edited and copied into the INCLUDE 
directory to become the new 
Name . INS 

files for the dynamically dimensioned 
UDFASTEX code. 

Development directory used for the new modal 
interpolation routines used in UDFASTEX. 
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MODE Contains all UDFASTEX Fortran codes, data, 

and procedure files for Option 3 shown on 
Table 5.1. 

PLOT Contains all UDFASTEX Fortran codes, data, 

and procedure files for Option 5 shown on 
T able 5.1. 

TAPE10 Contains all UDFASTEX Fortran codes, data, 

and procedure files for Option 4 shown on 
Table 5.1. 
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LISTING OF MAIN UDFASTEX PROCEDURE FILE : UD.COM 
DIRECTORY NAMES ARE THOSE USED ON THE UNIVERSITY OF DAYTON 
RESEARCH INSTITUTE VAX 11/780 


$ ! 

$ SET DEF USER DISK: [DRYDEN. RFT] 

$ ! 

$ ON CONTROL_Y THEN GOTO MAIN_MENU_NO_CLEAR 

$ ON ERROR THEN GOTO MA I N_MENU_N 0_CLEA R 

$ ON WARNING THEN GOTO MAIN MENU NO CLEAR 

$ ! - _ - 
$ TT :== WRITE SYS$OUTPUT 

$ CLEAR :== SET TERM / WIDTH = 80 

$ GOTOMAIN :== SET DEF [DRYDEN. RFT] 

$ GOTOINCLUDE : = = SET DEF [ DR YDEN. RFT. BIGGER. INCLUDE] 

$ GOTOMODE :== SET DEF [DRYDEN. RFT. REVISED. MODE] 

$ GOTOPLOT :== SET DEF [ DRYDEN. RFT. PLOT] 

$ G0T0TAPE1 0 : = = SET DEF [ DRYDEN. RFT . TAPE1 0 ] 

$ GOTOTEST :== SET DEF [DRYDEN. RFT. BIGGER. INCLUDE. TEST] 

$ ! 

$ CLEAR 
$ TT " " 

$ TT " " 

$ TT ” " 

$ TT " " 

$ TT ” " 

$ TT " " 

$ TT " " 

$ TT " " 

$ TT " 

$ tt " " 

$ tt " MASTER PROCEDURE ACTIVATED FOR RUNNING " 

$ TT " ALL REVISED FASTEX FLUTTER ANALYSIS CODES " 

$ TT " DEVELOPED BY THE UNIVERSITY OF DAYTON " 

$ TT " FOR NASA DRYDEN FLIGHT RESEARCH CENTER " 

$ TT ” UNDER GRANT NAG 2-377 " 

$ TT " ( LAST REVISED FEBRUARY 6, 1987 ) " 

$ TT " " 

$ TT " 

$ TT " " 

$ INQUIRE YN - 

" DO YOU WANT TO CONTINUE( Y/N) ? " 

$ IF ( YN .EQ. "Y" ) THEN GOTO MAIN_MENU 
$ GOTO QUIT 

$ j a#######*###########*####*#***##****#******#########*###*#**##**#*# 
$ MAIN_MENU : 

$ i ******************************************************************* 

$ CLEAR 

$ MAIN_MENU_NO_CLEAR: 

$ TT ” " 

$ TT " " 
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$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


TT w 

TT if 

TT ,T 

TT " OPTION 

TT " 1 

Tt if 2 

TT « 3 

'p-p tt 1( 

TT " 5 

TT " 6 

pp it 7 

pp II 

TT " 8 

TT tl 

tt it 9 

TT " 99 

pp II II 

INQUIRE OPT - 


UDFASTEX MAIN MENU 


DESCRIPTION 

Run UDFASTEX with internal mode interpolation” 
Run UDFASTEX with external mode interpolation" 
Compute aero box data or run plate spline" 

Run an alternate aero or eigen solution" 

Plot flutter V-g or V-f data" 

Change the problem size 
Compile and link an new 
for OPTION 1 use" 

Compile and link an new version of UDFASTEX" 
for OPTION 2 use" 

Create a print file of selected UDFASTEX output" 
Quit" 


capacity of UDFASTEX" 
version of UDFASTEX" 


ENTER OPTION" 


$ 

IF 

( 

OPT 

.EQ. 

II 1 II 

) THEN 

GOTO 

RUN INTERNAL 

$ 

IF 

( 

OPT 

. EQ. 

ti 2 « 

) THEN 

GOTO 

RUN EXTERNAL 

$ 

IF 

( 

OPT 

.EQ. 

it b it 

) THEN 

GOTO 

AERO SPLINE 

$ 

IF 

( 

OPT 

. EQ. 

II Ij II 

) THEN 

GOTO 

AERO EIGEN 

$ 

IF 

( 

OPT 

.EQ. 

It ^ ft 

) THEN 

GOTO 

PLOT - 

$ 

IF 

( 

OPT 

. EQ. 

"6" 

) THEN 

GOTO 

CHANGE SIZE 

$ 

IF 

( 

OPT 

.EQ. 

it j it 

) THEN 

GOTO 

COMPILE AND LINK1 

$ 

IF 

( 

OPT 

. EQ. 

II 8 H 

) THEN 

GOTO 

COMPILE AND LINK2 

$ 

IF 

( 

OPT 

. EQ . 

IT 9 II 

) THEN 

GOTO 

select "Print 

$ 

IF 

( 

OPT 

. EQ. 

II 9 9 TT 

) THEN 

GOTO 

QUIT 

$ 

TT 

Tt 


INVALID OPTION - 

- TRY 

AGAIN. . ." 


$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 


GOTO MAIN_MENU : 

! a#*###################*################*##########*##*######*#*#### 
a###****#########**#*##*##****##*#*###**********#**#**********#**#* 

RUN_INTERNAL : 

OPTION 1 Run UDFASTEX with internal mode interpolation 
********** ********************************************************* 

CLEAR 

pp ii it 

TT " OPTION 1 RUN FASTEX WITH INTERNAL MODE INTERPOLATION" 

TT " IS STARTING. . 

pp ii ii 

0PT1 : 

pp ii ii 

GOTOINCLUDE 
ON ERROR 


THEN GOTO MAIN MENU NO CLEAR 


@0PTI0N_1_REQUIREMENTS 

pp ii ii 

TT " DO YOU WANT TO" 

TT " 1) CHECK DIRECTORY" 

TT " 2) CONTINUE" 

TT " 3) GO BACK TO MAIN MENU" 

INQUIRE OPT - 

" ENTER OPTION" 

IF ( OPT .EQ. "1" ) THEN GOTO OPT1 1 
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$ IF ( OPT . EQ. "2" ) THEN GOTO OPT1 2 
$ IF ( OPT .EQ. "3" ) THEN GOTO OPT1 3 
$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN...*' 

$ GOTO OPT1 

$ OPT11: ! CHECK FOR FILES 

$ DIR .DAT 
$ DIR .EXE 
$ GOTO OPT1 

$ 0PT12: ! START PROCEDURES 

$ ON ERROR THEN GOTO MAIN_MENU_NO_CLEAR 

$ @START_OPTION_1 
$ GOTOMAIN 
$ GOTO OPT1 

$ OPT1 3 : ! GO BACK TO MAIN MENU 

$ GOTO MAIN_MENU 

$ j ******************************************************************* 
$ i ******************************************************************* 

$ RUN_EXTERNAL : 

$ ! OPTION 2 Run UDFASTEX with external mode interpolation 

$ j a##**###**#***##*##*####*#**########*#*###*######################## 

$ CLEAR 
$ TT " " 

$ TT " OPTION 2 RUN UDFASTEX WITH EXTERNAL MODE INTERPOLATION" 

$ TT " IS STARTING. . ." 

$ TT " " 

$ OPT2 : 

$ TT " " 

$ GOTOINCLUDE 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_CL E A R 

$ §OPTI ON_2_REQUI RE ME NTS 
$ TT " " 

$ TT " DO YOU WANT TO" 

$ TT " 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT21 

$ IF ( OPT .EQ. "2" ) THEN GOTO 0PT22 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT23 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT2 

$ 0PT21 : ! CHECK FOR FILES 

$ DIR .DAT 
$ DIR .EXE 
$ GOTO OPT2 

$ 0PT22 : ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_MENU_N 0_CLEA R 

$ @START_OPTION 2 
$ GOTOMAIN 
$ GOTO OPT2 

$ OPT2 3 : ! GO BACK TO MAIN MENU 
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$ GOTO MAIN_MENU 

$ t ******************************************************************* 
$ i ******************************************************************* 

$ ’ AERO_SPLINE: 

$ ! OPTION 3 Compute aero box data or run plate spline 

$ i a###*####*####**###*#########**#####*####*###*#**##*#######*####*## 

$ CLEAR 
$ TT " " 

$ TT " OPTION 3 COMPUTE AERO BOX DATA OR RUN PLATE SPLINE" 

$ TT " IS STARTING. . 

$ TT " " 

$ OPT3: 

$ TT " " 

$ GOTOMODE 

$ ON ERROR THEN GOTO MA IN_MENU_NO_CLEAR 

$ @OPTI ON_3_REQUI RE ME NTS 
$ TT " " 

$ TT " DO YOU WANT TO" 

$ TT " 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT31 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT32 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT33 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT3 

$ OPT31 : ! CHECK FOR FILES 

$ DIR .DAT 
$ DIR .EXE 
$ GOTO OPT3 

$ OPT32: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_C LE A R 

$ @START_OPTION_3 
$ GOTOMAIN 
$ GOTO OPT3 

$ OPT33: ! GO BACK TO MAIN MENU 

$ GOTO MA IN_MENU 

$ j ft****************************************************************** 
$ i ******************************************************************* 

$ AERO_EI GEN : 

$ ! OPTION 4 Run an alternate aero or eigen solution 

$ j ******************************************************************* 

$ CLEAR 
$ TT " " 

$ TT " OPTION 4 RUN AN ALTERNATE AERO OR EIGEN SOLUTION 
$ TT " IS STARTING. . ." 

$ TT " " 

$ OPT4: 

$ TT " " 

$ GOTOTAPE1 0 

$ ON ERROR THEN GOTO MAIN MENU NO_CLEAR 
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$ @0PTI0N_4_REQUIREMENTS 
$ TT " " 

$ TT " DO YOU WANT TO" 

$ TT " 1 ) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT41 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT42 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT43 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT4 

$ OPT41 : ! CHECK FOR FILES 

$ DIR .DAT 
$ DIR .EXE 
$ GOTO OPT4 

$ OPT42: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_MEN U_N 0_C LE A R 

$ @START_OPTION_4 
$ GOTO MAIN 
$ GOTO OPT4 

$ OPT43: ! GO BACK TO MAIN MENU 

$ GOTO MA IN_MENU 

$ ! ******************************************************************* 
$ j ******************************************************************* 

$ PLOT : 

$ ! OPTION 5 Plot flutter V-g or V-f data 

$ ; a####**#*#*###############***#*####*#**##########*#**###*#####*#### 

$ CLEAR 
$ TT " " 

$ TT » OPTION 5 PLOT FLUTTER V-G OR V-F DATA" 

$ TT " IS STARTING. . ." 

$ TT !t fr 
$ OPT5: 

$ TT f! n 
$ GOTOPLOT 

$ ON ERROR THEN GOTO MA I N_ME NU_NO_CLEAR 

$ @OPTION_5_REQUIREMENTS 
$ TT " ff 

$ TT " DO YOU WANT TO" 

$ TT " 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT51 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT52 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT53 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT5 

$ 0PT51 : ! CHECK FOR FILES 
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$ DIR FLUTTERPLOTFILE.DAT 
$ GOTO OPT5 

$ OPT52: ! START PROCEDURES 

$ ON ERROR THEN GOTO MAIN_MENU_NO_CLEAR 

$ @START_OPTION_5 
$ GOTOMAIN 
$ GOTO OPT5 

$ OPT53: ! GO BACK TO MAIN MENU 

$ GOTO MA IN_MENU 

$ ! a*##**#*#*###*###*######**###############*#*#*##*#*#**###*####**##* 

$ i a####*#*##########*##*####*###########*##*##***##*##**#*##*#*#*###* 

$ CHANGE_SIZE : 

$ ! OPTION 6 Change the problem size capacity of UDFASTEX 

$ i ******************************************************************* 

$ CLEAR 
$ XT " •' 

$ TT " OPTION 6 CHANGE THE PROBLEM SIZE CAPACITY OF UDFASTEX” 

$ TT " IS STARTING. . 

$ TT " " 

$ OPT6: 

$ TT " " 

$ GOTOTEST 

$ ON ERROR THEN GOTO MA IN_MENU_NO_CLE AR 

$ §OPTI ON_6_REQUI RE ME NTS 
$ XT ” " 

$ XX " DO YOU WANT TO” 

$ TT ” 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT61 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT62 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT63 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT6 

$ OPT 61 : ! CHECK FOR FILES 

$ DIR . $NS 
$ GOTO OPT6 

$ OPT62: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_C L E A R 

$ @START_OPTION_6 
$ GOTOMAIN 
$ GOTO OPT6 

$ OPT63: ! GO BACK TO MAIN MENU 

$ GOTO MAIN MENU 

$ ; ****** *T********************* ********** ******************* ********* 

$ ; ******************************************************************* 
$ COMP ILE_AND_LINK1 : 

$ ! OPTION 7 Compile and link an new version of UDFASTEX 
$ ! for OPTION 1 use 

$ i ******************************************************************* 

$ CLEAR 
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$ TT " " 

$ TT " OPTION 7 COMPILE AND LINK A NEW VERSION OF UDFASTEX1 " 

$ TT " IS STARTING. . 

$ TT " " 

$ OPT7: 

$ TT " " 

$ GOTOINCLUDE 

$ ON ERROR THEN GOTO MAIN MENU NO CLEAR 

$ @OPTION_7_REQUI RE ME NTS “ 

$ TT " " 

$ TT " DO YOU WANT TO" 

$ TT " 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 

$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT71 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT72 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT73 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT7 

$ OPT71 : ! CHECK FOR FILES 

$ DIR .FOR 
$ DIR .INS 
$ GOTO OPT7 

$ OPT72: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_C LE A R 

$ @START_OPTION_7 
$ GOTOMAIN 
$ GOTO OPT7 

$ OPT73: ! GO BACK TO MAIN MENU 

$ GOTO MAIN_MENU 

$ j a#*##############*##*##*###*#*#####*#*#####*#*####*##**########*### 

$ j ft##*###############*#########**######*##*#########*#####*#*#***#**# 

$ C0MPILE_AND_LINK2 : 

$ ! OPTION 8 Compile and link an new version of UDFASTEX 
$ ! for OPTION 2 use 

$ j a#*#####################*##*#****#**##*##*#*#*######**###*#######*# 

$ CLEAR 
$ TT " " 

$ TT " OPTION 8 COMPILE AND LINK A NEW VERSION OF UDFASTEX2" 

$ TT " IS STARTING. . ." 

$ TT " " 

$ OPT8: 

$ TT " " 

$ GOTOINCLUDE 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_CL E A R 

$ §OPTION_8_REQUIREMENTS 
$ TT " ,r 

$ TT " DO YOU WANT TO" 

$ TT " 1) CHECK DIRECTORY" 

$ TT " 2) CONTINUE" 

$ TT " 3) GO BACK TO MAIN MENU" 
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$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT . EQ. "1" ) THEN GOTO OPT81 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT82 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT 83 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT8 

$ OPT81 : ! CHECK FOR FILES 

$ DIR .FOR 
$ DIR .INS 
$ GOTO OPT8 

$ OPT82: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_MEN U_N 0_CLE A R 

$ @START_OPTION_8 
$ GOTO MAIN 
$ GOTO OPT8 

$ OPT 83 : ! GO BACK TO MAIN MENU 

$ GOTO MAIN_MENU 

$ I a#######*#******#***#*******##*#*#**#*#*#***#**#*#*###****#*#*#*##* 

$ SELECT_PRINT : 

$ ! OPTION 9 Create a print file of selected UDFASTEX output" 

$ J ****************************** ************************** *********** 

$ CLEAR 
$ TT " " 

$ TT " OPTION 9 CREATE A PRINT FILE OF SELECTED UDFASTEX OUTPUT" 

$ TT " IS STARTING. . ." 

$ TT " " 

$ OPT9 : 

$ TT " " 

$ GOTOINCLUDE 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_CL E A R 

$ §OPTION 9 REQUIREMENTS 


$ 

TT 

TT 

TT 


$ 

TT 

TT 

DO YOU WANT TO" 

$ 

TT 

TT 

1 ) 

CHECK DIRECTORY" 

$ 

TT 

TT 

2) 

CONTINUE" 

$ 

TT 

TT 

3) 

GO BACK TO MAIN MENU" 


$ INQUIRE OPT - 

" ENTER OPTION" 

$ IF ( OPT .EQ. "1" ) THEN GOTO OPT91 

$ IF ( OPT .EQ. "2" ) THEN GOTO OPT92 

$ IF ( OPT .EQ. "3" ) THEN GOTO OPT93 

$ ! INVALID ENTRY 

$ TT " INVALID OPTION - TRY AGAIN..." 

$ GOTO OPT9 

$ OPT91 : ! CHECK FOR FILES 

$ DIR UDFAS*.LIS 
$ GOTO OPT9 

$ OPT92: ! START PROCEDURES 

$ ON ERROR THEN GOTO MA I N_ME N U_N 0_C LEAR 

$ @START_OPTION_9 
$ GOTOMAIN 
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$ GOTO OPT9 

$ OPT93 : ! GO BACK TO MAIN MENU 

$ GOTO MA IN_MENU 

$ ; ******************************************************************* 

$ i ******************************************************************* 

$ * QUIT: 

$ ! OPTION 99 Quit 

$ ! ******************************************************************* 
$ GOTOMAIN 
$ EXIT 
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APPENDIX C 


PLATE SPLINE THEORY 



BIVARIATE DATA INTERPOLATION 


The following is based on the theory presented in Reference 12. 
1. Usage - 

a. Calculation of modal data in aerodynamics grid. 

b. Calculation of aerodynamic matrices at intermediate 
values of reduced frequency and Mach number. 




C-l 



3. Equations - 

Assume the plate spline function in the form 

N 

(1) W(x,y) = a Q + a^ + a 2 y + ^ K i (x,y)P i 
where 


(2) (x,y) = r i 2 lnr i 2 

(3) r. 2 = (x - x. ) 2 + (y - y ) 2 

^ X 1 


and a^ , a^, a 2 , are constants to be determined by 

the conditions 


N 

(4) l P = 0 
i=l 1 


N 

(5) l x P = 0 
i=l 1 1 


N 

(6) l y P = 0 
i=l 1 1 

N 

(7) Wj = a Q + a x Xj + a^ + t X K^Pj. ,j=l,...,N 


i=j 


where 


(8) K. . = 
13 


K i < x j 'Vj ) 


Note K. . = K . . . 

T -I -IT 
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Write (.4) - (.7) in the matrix form of N+3 equations in 
N+3 unknowns: 


C9) 


I A] {c} = ib} 


where the spline matrix is 


(10) [A] = 


0 0 0 

0 0 0 

a o o 


i 


1 -X, 


1 X. 


N J N 


A 1 

*1 


K 


21 


K 


12 

0 


N 


N 


K 


IN 


K 


2N 


K N1 K N2 • • * ° 


N+3 by N+3 


Cll) 


(b> = < 


0 

0 

0 

W ] 

w. 


w. 


V N J 


N+3 by 1 


and the desired solution is the set of spline coefficients 




( 12 ) 


(c) = { Pj. 




V P N / 


N+3 by 1 
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The required derivatives of W are 


(13) a l + J 1 K xi (x 'l' )P i 


N 


where 


(14) K xi (x,y) = 2 (x - x^ ) (1 + lnr^ ) 


and 


(15) 

where 


3W 


N 


-s— = a_ + £ K . (x,y)P. 

9y 2 i^ yi ' * 1 


(16) Kyi ( x / y ) = 2 (y - y^ (1 + lnri ) 


4. Computational Procedure 


STEP 


COMMENTS 


1. Input N 


N = Number of given data points 


X 1 ' ‘ * ' ' X N 
y l' * * * ,y N 


Xw...,x = x coordinates of 

given data 

y l'***' y N = y coordinates of 
1 given data 

W. , . . . , W = Given functional 
1 N values 


2. Input M 
X 1 ' * ' * ' x m 
y l' * * * ' y M 


M = Number of points at which 

W, W , W are needed 
x y 

x.,...,x = x coordinates of 

needed points 

y l''*‘' y M = y coor, 3inates of 
1 needed points 


3. Compute by 

Equation (8) 


= Spline weights. Compute 
1 - 5 Kij = 0 separately, 

In terms become large 
otherwise. 
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I 


STEP 

4. Set up [A] matrix by 
Equation (10). 


5. Solve Equations (.3) 
for Ic} 


6. Use {c} in Equation (.1) 

to compute 

7. Use (c}_in Equation (.13) 

to get W x;j 

8. Use {c}_in Equation (15) 

to get 


COMMENTS 


Note IAJ is symmetric 


If Gauss elimination is used, 
note that diagonals are zero. 
Use symmetry of I A] in solution 
process if possible. 


Wj = desired functional values 


W . = x derivatives desired 

xd 


W . = y derivatives desired 
YJ y 
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5. Extensions 


Improved accuracy may he obtained by coordinate transforma- 
tions on the Cx.,yji, U_.,y..) grids. The following steps 
are proposed: 

1. Translation of data to centroid 

(17) V " * *0 

U8) yj ' = yj - y 0 

where 

N 

U9) X “ = B 

(20) y 0 = N //j 

2. Rotation of data to principal axes 


(21) 

x. " 

= x . 1 cos8 

+ y^' 

sin© 

3 

3 

3 


(22) 

V 

= x . ' sin0 
3 

♦ V 

COS0 

where 


1 j 

' -21 , 

1 fl 

(23) 

8 = 

1/2 tan _i ( 

x 1 

y 

I , - 

i 1 



V 

■ x' 

y * 


(Note-use ATAN2(-2l x , y , , (I x « - I y ,))) 
N 2 

(24) I , = l (y,’) 
x j=l 3 
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(25) 


(26) 



N 

l Cx'l 
j=l 3 


x'y ' 


N 

l u.')(y,') 

j=i 3 3 


3. Scaling of data 


(27) 

x. m = 
3 

x "/x" 
x j / max 

(28) 

y = 

J j 

y ."/y" 

2 3 ' 2 max 

where 




(29) 

ti _ 

max 

max 

j=l,N 

| x . " | 
3 

(30) 

y " = 

2 max 

max 

j=l,N 

1 y j " 1 


Apply Equations (17) - (30) to (x^y^), (x^y^) then go to 

Step 3 of computational procedure. 


(x. m 

»Yj ) 

in place 

of 

< 3 V Y j> 

(Xj'" 

#Yj ) 

in place 

of 

(X j ' y j* 
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